Hybrid Levenberg–Marquardt and weak-constraint ensemble Kalman smoother method
نویسندگان
چکیده
The ensemble Kalman smoother (EnKS) is used as a linear least-squares solver in the Gauss–Newton method for the large nonlinear least-squares system in incremental 4DVAR. The ensemble approach is naturally parallel over the ensemble members and no tangent or adjoint operators are needed. Furthermore, adding a regularization term results in replacing the Gauss–Newton method, which may diverge, by the Levenberg–Marquardt method, which is known to be convergent. The regularization is implemented efficiently as an additional observation in the EnKS. The method is illustrated on the Lorenz 63 model and a two-level quasigeostrophic model.
منابع مشابه
A hybrid (variational/Kalman) ensemble smoother for the estimation of nonlinear high-dimensional discretizations of PDE systems
Two classes of state estimation schemes, variational (4DVar) and ensemble Kalman (EnKF), have been developed and used extensively by the weather forecasting community as tractable alternatives to the standard matrix-based Kalman update equations for the estimation of high-dimensional nonlinear systems with possibly nongaussian PDFs. Variational schemes iteratively minimize a finite-horizon cost...
متن کاملControlling balance in an ensemble Kalman filter
We present a method to control unbalanced fast dynamics in an ensemble Kalman filter by introducing a weak constraint on the imbalance in a spatially sparse observational network. We show that the balance constraint produces significantly more balanced analyses than ensemble Kalman filters without balance constraints and than filters implementing incremental analysis updates (IAU). Furthermore,...
متن کاملA Hybrid Monte-Carlo Sampling Smoother for Four Dimensional Data Assimilation
This paper constructs an ensemble-based sampling smoother for fourdimensional data assimilation using a Hybrid/Hamiltonian Monte-Carlo approach. The smoother samples efficiently from the posterior probability density of the solution at the initial time. Unlike the well-known ensemble Kalman smoother, which is optimal only in the linear Gaussian case, the proposed methodology naturally accommoda...
متن کاملA Carrier Estimation Method Based on MLE and KF for Weak GNSS Signals
Maximum likelihood estimation (MLE) has been researched for some acquisition and tracking applications of global navigation satellite system (GNSS) receivers and shows high performance. However, all current methods are derived and operated based on the sampling data, which results in a large computation burden. This paper proposes a low-complexity MLE carrier tracking loop for weak GNSS signals...
متن کاملThe Iterative Ensemble Kalman Smoother: the Best of Both Worlds?
Data assimilation seeks a mathematically optimal compromise between outcomes of a numerical model that simulates a physical system and observations of that system. It has been successfully used for twenty years in operational meteorology to perform the best forecast, and is now being used or tested in many geoscience fields. Two main classes of methods have taken the lead. Firstly, 4D-Var is a ...
متن کامل